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1 Introduction 

Ordinary matter is composed of electrons (negatively charged) and nuclei (positively 
charged) interacting via electromagnetic forces. The electric potential between two 
such particles of charges Q and Q' located at r and r' in R3 is ee7|r-r'|. The 
Coulomb potential poses two difficulties: (1) The local singularity and (2) its long 
range. One has to understand why the local singularity does not cause instabilities 
and why the long range does not have strong macroscopic effects. One of the first 
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major triumphs of the theory of quantum mechanics is the explanation it gives of 
the stability of the hydrogen atom (and the complete description of its spectrum) 
and of other microscopic quantum Coulomb systems. It, surprisingly, took nearly 
forty years before the question of stability of everyday macroscopic objects was 
even raised. The rigorous answer to the question came shortly thereafter in what 
came to be known as the Theorem on Stability of Matter proved first by Dyson 
and Lenard and later, in a more simple and transparent way by Lieb and Thirring. 
Since these seminal works, the proof of stability of matter has been extended to 
several different settings, including relativistic systems and systems in the presence 
of dynamic electromagnetic fields. 

We will, in particular, discuss the importance of particle statistics, i.e., whether 
the particles are bosons or fermions. Using Bogolubov's theory for Bose gases 
Dyson concluded that charged bosonic particles would be macroscopically unsta- 
ble. That Bogolubov's theory gives the correct description of charged Bose systems 
was proved in a series of papers by the author in collaboration with E.H. Lieb. 

Having proved stability of matter the next question is whether one can establish 
the thermodynamics of charged systems, i.e., the existence of the thermodynamic 
limit. This was originally achieved by Lieb and Lebowitz. We will describe a new 
approach to the existence of the thermodynamic limit, which applies to many dif- 
ferent quantum Coulomb systems, possibly in the presence of an underlying lattice 
structure. This generalizes earlier work of Fefferman. 

An important theme in these notes is the use of functional inequalities, among 
which a prominent role is played by the Lieb-Thirring inequality. 

Another important theme will be how to control the screening of the Coulomb 
potential. 

The notes are organized as follows. In Section |2] we derive the classical Hamil- 
tonian for charged particles interacting with electromagnetic fields. In Section[3]we 
discuss different ways of quantizing the system. We may quantize the particles and 
leave the fields unquantized or quantize both the particles and the fields. Moreover, 
we may consider the particles relativistically or classically. This will give a variety 
of different models. In Section |4] we discuss stability both in the sense of stability 
of individual atoms and in the sense of stability of macroscopic matter Finally, in 
Section|5]we discuss situations where stability fails, in particular the case of bosonic 
matter. 



2 Classical "point" charges 

We consider classical particles with charges Qit.-tQn G K situated at points 
ri , . . . , r^i G E^. Dimension 3 is of course the physical space dimension, but a natu- 
ral questions would be whether the discussion of charged system could be general- 
ized to other dimensions. It is however not entirely clear what the correct physical 
questions would be and we will there restrict the discussion to dimension 3. 
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We would ultimately like to consider the particles as point charges. Unfortu- 
nately this will however lead to divergencies. In order to avoid this we will ini- 
tially assume that each particle is given by a charge distribution QiXR{r — r,) where 
Xr{'') =R^^X{>'/R) X G Co(IR^) (a continuous compactly supported function) 
with/Rs;); = 1. 

We will start the discussion of this system of charged particles from the La- 
grangian of the particles and the electromagnetic field. It i|3 (using Gaussian units) 

L«(r,,A,V) = f {T*{rj)-QM*XR{rj)-QjV*XR{r,i)) 

where 

• Vj denotes the velocity of particle j. 

• (v) is the Legendre transform of the kinetic energy Tj{p) as a function of mo- 
mentum p. We assume that the kinetic energy functions Tj are convex functions 
on K^. For a non-relativistic particle of mass nij we have Tj{p) = j^P'^ ™d thus 

T^(y) = jnijv^ and for a relativistic particle it is Tj{p) = \J + ntj — mj and 

hence TJ (v) =mj — mjy/ 1 — v^. We are using units in which the speed of light is 
1 (a convention we will use throughout these notes). 

• A is the vector potential and the magnetic field is B = V x A 

• y is the electric potential. The electric field is 

E = -<9,A-vy. 

In order to go to a Hamiltonian description we will choose Coulomb gauge or more 
precisely require that 

V • a,A = 0, 

i.e., we assume that V - A is independent of time. We will also for simplicity assume 
that A decays sufficiently fast that we are allowed to ignore boundary terms when 
integrating by parts. 

With the Coulomb gauge condition we see that 

fix = -d,A 

is the divergence free part of the electric field. The total electric field is £ = Ej_ — W 
and we have 

||£|2 = j\E^\^ + j\VV\\ 

' Strictly speaking, even if we use a relativistic kinetic energy, this Lagrangian is not relativistically 
invariant. The reason is that we consider the particles as rigid bodies, which do not Lorenz contract 
as they move. We will here ignore this additional complication. The Lagrangian in the form given 
here is that of the Abraham model of charged particles 1 321 . 
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The electric potential V is not a dynamic variable in the sense that that dtV = V 
does not occur in Lr. The equation for V is 

where we have used the Coulomb gauge condition. We recognize the equation for 
V as Gauss' law 

N 

4?r £ QjXR{r-rj) = -AV = -V • {dtA + VV) = V • £. 
;=i 

In the Hamiltonian formaUsm this is a constraint equation. The solution is 

N 



V{r) = '£Qj [xR{r)*\r-rj\-'dr. 



The canonical variable dual to rj is 



Pj = WrjLR = WT* (o) - QjA * XRin). 
The canonical variable dual to A is 



^-—d,A = -{4n)-'Ej 



5d,A An 

due to the Coulomb gauge condition. We then find the Hamiltonian function (where 
we assume that Tj is convex such that the double Legendre transform of TJ* = Tj) 

N 1 , 

HR{rj,pj,A,E^) = Y^pji-j-— J dtA-E^-LR{rj,A,V) 

N N 

= L Tjipj + QjA * XR{rj) ) + £ QjV * XR{rj) 

+^/(|ExP + |VxAP)-±/|vyp 

= L Tj [pj + QjA *XR(rj))+^Y. Qj^ * XR{rj) 

7=1 ^ ^ ^7=1 

+ 4/(I^.P + |VxAP) 
ZTj{pj + QjA*XR{rj)) 



7=1 
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If we subtract the (divergent) self-energy 

Pi2 J J |r-r'| 7^1 2 I'"- '"'I 



and assume that A is continuous at rj, we see that as R tends to zero Hr converges 
pointwise to the Hamiltonian 

i/(o-,p,-,A,£j.) = £7}(p,+eAo))+ I 

+ + (1) 

Here r^- and pj are dual canonical variables as are the field variables A and 
— {4n)~^E±. Unfortunately, this Hamiltonian is only formal and suffers from sin- 
gularities (the field A that solves Hamilton's equations will be singular at rj) leading 
to severe difficulties in describing the motion of classical point charges. 
If external fields Vex and Aex are present, the energy is 



I {tj{pj + Qj{A +Aex)(o)) + QMrj)) + ^^Y._^ML. 

IJ\\E^\^ + \B\'). 



j=l " ' l<i<j<N I'"' 



3 Charged quantum gases 

We now discuss the quantization of the Hamiltonian of charged point particles. We 
emphasize that it has not been possible to define a fully relativistically invariant and 
causal theory of quantum electrodynamics (QED) and all the models we describe 
here are at best approximations to such a theory (if it exists). All models discussed 
here are mathematically well defined (except when otherwise stated expUcitly). 
There are several levels of quantization that may be considered. 

• We can leave the fields A and E± classical and quantize the particles, i.e., de- 
scribe them by a square integrable wave function i/A(ri , . . . , r^). 

• We may quantize the particles and the fields, i.e., describe the particles in terms 
of a wave function and turn the fields into operator valued functions A and E±. 
This would require introducing some cut-off regularization in the fields. 
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• We may second quantize the particles, i.e., let also y/^ be an operator valued func- 
tion. This procedure is necessary if we consider relativistic particles described by 
the Dirac operator. 



3.1 Quantized particles and classical fields 

The variables of the system are the 3-dimensional vector fields A,Ej^ (assumed to 
satisfy appropriate regularity and decay properties, at least implying that and 
B = V X A are square integrable) and the wave function which is a normalized 
function in (g)^2[L^(i2)]^>, where ^2 c (say an open set) and Vj is a positive 
integer counting the number of internal degrees of freedom of particle j, (e.g. a 
particle of spin s would correspond to = 2i + 1). We shall write 

Y=Y{n,si,...,rN,SN), rj e Q, sj = l,...,Vj. 

The energy of the system is given by 

^N{W,A,E) = {xif,HN{A,E^)xi/), (2) 

where (v^, 0) refers to the inner product of S 0^ JL'^(X2)]^j and Hi^ is the 
(unbounded) operator (depending on A and E±) 

HN{A,E^)=l^(Tj{-iVj + Qj{A+A,,){rj))+QjV,,{rj))+ £ 

+ ^J\\EA^ + \B\^ (3) 

We will throughout be using units in which the reduced Planck constant h= I. 
The last integral above acts as a (A and E± dependent) scalar in the Hilbert space. 
The Hamiltonian H{A,E±) depends also on the exterior fields Vex and Aex, but we 
suppress this in the notation as these fields usually remain fixed. In fact, we will 
mostly, and unless otherwise exphcitiy stated, assume that the exterior fields vanish, 
i.e.. Vex = and Aex = 0. 

The expectation value S'n{y,A,E±) is not defined for all ^ in the Hilbert space 

.7=1 

We will here avoid the discussion of domains of self-adjointness of the operator 
H{A,Ej^). We will instead restrict attention to ^j/ in the the subspace of smooth 
functions with compact support, i.e, [Cq (f2^)]^i- '^'^' C (gj^^j [L^{Q)YJ. 

One of the main issues we will discuss in these notes is the question of stability, 
i.e., whether <^(y/^, A, is bounded from below independently of xj/ (normalized). 
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A , and E±. If such a lower bound holds the operator H{A,E±) has a self-adjoint 
Friedrichs extension and we are actually making claims about this extension. From 
our point of view the only complication due to considering the restriction to Cq is 
that a possible ground state (a state achieving the lowest possible energy) is most 
likely not represented by an element in Cq, but only by an element in the Friedrichs 
extended domain. We shall, however, not be concerned with the actual ground states, 
but only the energy, so we ignore this issue. 

Since the three coordinates of —i'Vj + QjA{rj) correspond to in general non- 
commuting operators, we must discuss the meaning of 7)( — iV j + QjA{rj)y We 
will, in fact, only consider examples where the functions Tj{p) can be written in 
terms of (possibly matrix-valued) polynomial expressions of p in such a way that 
the meaning of 7} (at least on a suitable domain) will be clear. The examples we will 
consider are 

• Non-relativistic kinetic energy operators, where Tj{p) — {Inij)^^ p^, i.e., the op- 
erator is 

Tj{- /V, + QjA{rj)) = {2mj)-'(-Nj + QjA{n)f- (4) 

We will refer to particles with this kinetic energy as non-relativistic parti- 
cles. This is the kinetic energy used when treating non-relativistic atoms and 
molecules or ordinary matter 

• Relativistic kinetic energy operators, where Tj{p) = {p^ +mjy/^ ^"^y- the 
operator is 

Tj{-iVj + QjAirj)) ^ ((-/V, + e,A(r_,-))2+m2)i/2_„,^, (5) 

The square root of an operator is here defined in the spectral theoretic sensfi 
We will refer to particles with this kinetic energy as relativistic (or sometimes 
pseudo-relativistic) particles. Both relativistic and non-relativistic particles may 
have internal degrees of freedom corresponding to V/, being greater than one. 

• The non-relativistic and relativistic Pauli-operators. These are operators acting 
on two-component vector valued functions given by inserting the operator 

aj-i-Nj + QjA{rj)) 

into the kinetic energy functions above. Here CT = {<j\(7^ ,<J^) is the vector of 
2x2 Pauli matrices 

'■K-)' "^K"o> "'-{l-^)- 

(The subscript /' on (7 above indicates that it acts on the internal degrees of free- 
dom of particle j.) Thus in this case Vj~2. The resulting kinetic energy operators 
are 

Tj{-N, + Q,A{rj)) = (2m,)-'(<T,/- (-'V; + e;^(0)))' (6) 



^ The operator inside the square root is defined as a self-adjoint operator by Friedrichs extending 
it from the domain of smooth functions with compact support. 
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for non-relativistic PauU particles and 

1 /2 

Tj{-iyj + QjA{rj)) = (^{aj-{-Nj + QjAirj))y + mj) -mj (7) 

for relativistic Pauli particles. For the Pauli operator we have the Lichnerowicz' 
formula 

(ct . . (_,v,- + QjA{rj))^ = i-Nj + QjA{rj)f + QjGj ■ B{rj) (8) 

and we see that the Pauh operator includes the couphng of the particle spin to the 
magnetic field. 
• We could also consider the 4 x 4 Dirac operator 

Tj{p) = a-p + mjl5, 

.i.e., 

Tj{~Nj + QjAin)) - Uj ■ i-Nj + QjA{rj))+mjpj 
where a and j3 are standard 4x4 Dirac matrices, e.g.. 



a 



Oa\ 10 

<j r ^ [o-i 



using a 2 X 2-block notation. Thus in this case Vj = 4. In contrast to the other 
types of operators the Dirac operator, however, is not positive, in fact, not 
bounded below and we will therefore not be able to treat it unless we second 
quantize the particle fields. We will discuss this briefly below. A different ap- 
proach to deal with the unboundedness from below of the Dirac operator is to 
restrict to the subspace of L^(R^)^ which corresponds to the positive spectral 
subspace of the Dirac operator. This approach is caUed the no-pair theory, but we 
will not discuss it further here. 



3.2 Statistics of identical particles 

Until now all particles have been considered as distinguishable, but if we have iden- 
tical particles the issue of particle statistics plays an important role. 

The A^-particle space for A^-identical particles moving in X2 C and with V 
internal degrees of freedom is JifN — [L'^{^)V- On JlfN we define the orthogonal 
projections 

{P^\l/){ri,si,...,rN,SN) = -j^ £ (±l)''V^('-a-i(l)'*<T-i(l)'---''"<T-i(A')'*<T-i(A'))- 
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They project onto the symmetric (+) or antisymmetric (— ) subspaces. We denote 
these subspaces = P^Jf^. We will also use the notation — /\^[L^(i2)]*'. 

In case of identical particles, i.e., if the operators Tj and the charges Qj ~ Q 
are the same for all y = 1 , . . . , A^, the Hamiltonian Hn in (O maps the subspaces 
to themselves and it makes sense to restrict to these subspaces. We will write 

^±(VA,A,£x) (9) 

to emphasize that we restrict to e -y^j^- In the symmetric case (+) we say that 
we have a system of bosons in the antisymmetric case (— ) we say that we have a 
system of fermions. 

It is of course also possible to have mixtures of several species of identical parti- 
cles being fermions or bosons or even several species of identical particles together 
with a number of distinguishable particles. We leave it to the reader to work out the 
structure of the underlying Hilbert space and the Hamiltonian in the general case. 
We will look at specific examples later. 



3.3 Grand canonical picture 

It is often useful to consider a situation where the particle number is not specified 
at the outset, but where we would instead ask what the optimal particle number 
is in a given situation, e.g., what number of particles minimizes the energy. This 
picture is referred to as the grand canonical picture. The optimal particle number 
may if necessary be adjusted by adding a term fi times the particle number to the 
Hamiltonian. Such a parameter ji is called a chemical potential. We note that this 
is not the same as adding a constant to the exterior electric potential Vext as such a 
constant will multiply the total charge of the system. 

In order to treat variable particle number we define the bosonic or fermionic Fock 
spaces 

oo 

with the convention that J^f^ = C. The element 1 e C = Jff^ is referred to as the 
vacuum vector and will be denoted by |0) . For a normalized vector xj/ G we may 
write Xj/ = ©/^=o Wn, where y/N G with Y!n=o II VwIP = 1- We say that such a 
vector represents a grand canonical state and we define the grand canonical energy 
(with chemical potential included) 

.?±(M,V^,A,£x)- Y.'^NiVNA,E^) + ^iN\\wNf. (10) 

As before this energy is not defined for all but we restrict to Xj/ corresponding 
to finitely many particles, i.e., \j/ = 0^=o ¥n for some finite integer M and where 
each y/fj is in Cq. 
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Again it is possible to consider several species of identical particles in the grand 
canonical picture. We leave it to the reader to write down the Hilbert space and the 
energy (see also below). 



3.4 Second quantization and quantization of fields 

We shall here give a brief introduction to second quantization and discuss how to 
quantize particle fields and the electromagnetic fields. 

For f GL?{Qy we define the annihilation operator a (/) : Jif^ ^t^- 1 f or AT = 
l,...by 

(«(/)V')('"i>-si,---,'"iv-i,-SiV-i) = \ f{rN,SN)v{n,si,..-,rN,SN)drN. 

The adjoint of this operator is a* (/) : Jifk-i given by 

{a'{f)\^){ri,s„,...,rN,SN) = \/Nf{rN,SN)\^{n,si,...,rN-i,SN-i). 

We define the bosonic and fermionic annihilation operators a±{f ) : M'^ '^n-\ 
as the restriction of a{f) to the respective subspaces, i.e., a±(/) = a{f)^j^±. The 
adjoints are a^f) : ^ given by aHf) = P^a*{f\^±_^ . 

We may extend a± (/) and a±(/) to operators on the subspace of the Fock spaces 
corresponding to finite particle numbers, i.e., spanU^^Q ®n=(.)'^^n^- They can- 
not be extended as bounded operators on the full Fock spaces. 

The extended operators satisfy the famous conmiutation (+) and anti-commuta- 
tion (— ) relations 

[a±(/),ai(^)]± = {f,g)L2(nyI 

where [A , B]± = AB =F BA and / is the identity of Fock space (or rather its restriction 
to the subspace corresponding to finite particle numbers. 

If {fj} is an orthonormal basis in L^(i2)^ we define the operator valued distri- 
butions 

<t>±{r,s) = £/7('-,*)a±(/7), <Pl{r,s) = £ /X'-,*)«±(/j), 
allowing us to write 

«±(/) = L / 7M0±(n*Vr, al{f) = £ / f{r,s)r±{r,s)dr. 

s=\-' s=l-' 

for /eL^ (12) ^.Formally we have 
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We refer to and (j)* as field operators. 

Using tliese field operators we may write the grand canonical Hamiltonian 
®n=o{Hn + M^)> corresponding to identical fermions or bosons, formally as 

°° ^ /■ r 1 

^{Hn{A,E^)+HN)= £ r±ir,s) r^(-jV, + eA(r))+M5,y (|)±(r,/)Jr 

N=0 i,s'=K 

where 7"^^', s,/ = 1 , . . . , v refer to the matrix components of the kinetic energy op- 
erator. It is left as an exercise to the reader to check that the ordering of (j)* and (j> 
exactly gives the correct Hamiltonian with no-self interactions. 

In this formalism it is easy also to write down the grand canonical operator cor- 
responding to K different species of either fermions or bosons. For j = 1,...,K let 
Tj, Qj, and Vj represent the kinetic energy function, the charge, and the internal 
degrees of freedom of species j which is either a fermion or a boson. The relevant 
Hilbert space is = (S)f=i ^j{[L^{n)Vj) where .^j is the Fock space for species 
j. Denoting the field operators for species j by and (^^ the corresponding grand 
canonical Hamiltonian (with chemical potential fJ, included) is 

H{n,A,E^) =t L I \Tf{-iVr + QjA{r)) + n5J (^j{r,s')dr 

+i I I i // <^*{r,s)(l>J{r',s')^^j{r',s'Mr,s)drdr' 

^ ij=ls=\s'=l-' ■' l'^ I 

The subspace on which H [ji, A, E±) is defined is the space corresponding to finitely 
many particles and where the restriction to each particle sector is a smooth func- 
tion with compact support. Although it is hopefully clear what this means it is 
rather complicated to write it down explicitly. For the convenience of the reader 
we will nevertheless do this now. The subspace of Jff corresponding to Nj particles 
of species j, ,j = l,... ,K is (g)f=i Pj (0^^ [L'^{Q)Yj), where Pj refers to the rel- 
evant projection corresponding to the statistics of species j. We may consider this 

space a subspace of [L^(i2'^'+' '+'^'^)]^i "'^i^ . The subspace of smooth functions 

with compact support is [Cq (i2'^i+ '"^k . Thus the subspace on which we 
define //(/i,A,£j^) is 
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Ml Mk ( K Nj \ 

Mi,...Mj=QNi=Q Nk=Q \i=l / 

(11) 

The energy of the system with particles being in a state represented by i/A € ^ is 
denoted 

^ifi,xi/,A,E^) = {xi/,H{^i,A,EM- (12) 

It is easy to check that this agrees with the definition (fTOl i in the case of only one 
species. 



3.5 Quantization of the electromagnetic field 



We will briefly discuss how to quantize the electromagnetic field. We will remain in 
Coulomb gauge and quantize such that V A — Q. This is most conveniently done in 
momentum space. 

For^G M? choose ei (A:), 62 (^) such that ei (fc),e2(^),^ form an orthonormal 
basis. 61,62 cannot be chosen continuously, but this will not cause problems for what 
we want to say. 

Let (r, A ), A = 1 , 2 be a bosonic field operator with two internal degrees of free- 
dom. They are field operators for the light quanta, i.e., photons. Define the Fourier 
transformed operators (Of course they are also simply bosonic field operators) 

0(A:,A) = {2n)-^^^J e-''"'^{r,X)dr, 0*(/t,A) = {In)-^'^ j e''^''^* {r,X)dr. 

We define the quantized magnetic vector potential as the operator valued distribution 



Air) = (2k) 



-3/2 



X=\.2 



6^ [k) {e'^'ip {k,X)+ 6-'^'^ {k,X))dk (13) 



and the transversal electric field as 

-1/2 / . / 1^1 



E^ir) = i{2n)- 



1 = 1,2 



6^ik)i6'"-'^ik,X) 



-ikr 



(^{k,X))dk. (14) 



We then find the commutator between the conjugate variables 



Mr),-^E^,j{r'r) 



where P(r, r') is the 3 x 3-matrix valued integral kernel of the projection in L^{M?)^ 
projecting onto divergence free vector fields. 

A straightforward (formal) calculation gives for the field energy 
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i- f \EAr)\' + \VxAir)\^dr^\Y^f mrikAmA) 

on Jr3 2^jRi 

+${k,X)^*{k,l))dk. 

This expression however is infinite and we must normal order it to get a well-defined 
operator: 

y f \k\^*{k,l)^{k,X)dk. 
This is the field energy operator of the electromagnetic field on the Fock space 



3.6 Non-relativistic QED 

We may now write down the Hamiltonian of non-relativistic QED, i.e., of the quan- 
tized electromagnetic field coupled to quantized non-relativistic particles. The par- 
ticles will be described by the non-relativistic kinetic energies (01 or (|6]l, but since 
A is now an operator valued distribution, these operators will not make sense un- 
less we again introduce the extended charge distribution of the particles. The grand 
canonical non-relativistic QED Hamiltonian for K species of identical particles is 
then (ignoring for simplicity the chemical potential) 

^=E E U*{r,s)Tf{~N, + Q,A*XR{rmMs')dr 

+ E / \k\^*{k,X)^{k,X)dk. 
This operator is defined on the Hilbert space 

K 



The operators ^/ are field operators for the particles and ^ is the field operator for 
the photons. The energy may be calculated in a state represented by a 'f' in the 
subspace of the Hilbert space consisting of functions of finitely many particles 
and photons (we will not write this explicitly this time). The energy is denoted 



(15) 
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As written now the model depends on the regularization parameter R. The limit as 
R tends to is not well understood and will require at least to renormalize the bare 
mass and charges of the particles. 



3.7 Relativistic QED Hamiltonian 

As akeady emphasized a Hamiltonian (or for that matter any non-perturbative) for- 
mulation of QED is non-existent. Here we simply write down the formal expression 
for the Hamiltonian for the electron-positron field (with charge e) interacting with 
the electromagnetic field: 



Here <pe refers to the fermionic field operator for the electron-positron field and ^ 
is the bosonic field operator for the photon field. The operator A is given by (fTSl l. 
Note that we have not distinguished between electrons and positrons, but that the 
operator is written in a charge conjugation invariant way as the density is written as 

the commutator \Y^a=\ [0e(a,''),0e*(^j'')]+- 

The operator Hqed is ill-defined unless regularizations are introduced and even 
in this case it is very difficult to analyze. The no-photon situation was studied in the 
mean-field approximation in |15J . 



In the previous section we discussed how to define the energy of states of charged 
quantum gases in different models. 

We have introduced the fixed particle number (or canonical) energy (oA?(y/,A,Ex) 
in (|2|i (or the bosonic or fermionic analogs in (|9]l) or the grand canonical en- 
ergy ff (jti; </,A,£x) in (fT2T l. We also defined the non-relativistic QED energy 
<5NRQED(f) inO- 

We will say that a system is stable of the first kind or canonically stable if the 
energy <o]v(v,A,£x) is bounded below independently of A, Ei_, and normalized 
i/. In this case we will call the infimum of S'N{y/,A,E) the ground state energy 
regardless of whether an actual minimizer (a ground state) exists or not. Thus the 
canonical ground state energy of the system is 
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N 

£w(i2)=inf{jW(r,A,E) I i/e (0L^(r2)^.')nCo~(f2^)^'-^", ||r|| - 1, 

.7=1 

A,Ex eCo(K^;M^)}. 

Note that we are restricting the particles to be in the set £2 whereas A and E are 
unrestricted vector fields in M?. It is immediate to see that we might take Ej_ = in 
the infimum, this will however not be the case for quantized fields below. 

The ground state energy of course depends on the types of particles in the system. 
We are suppressing this dependence in order not to overburden the notation. 

The ground state is the state of the system at absolute zero temperature. It is of 
course also of interest to study quantum gases at positive temperature corresponding 
to minimizing the free energy we shall however not do this here. 

We could also have chosen to consider the purely static Coulomb potential and 
set A = 0, but as we shall see the inclusion of A does not really change the treatment 
in the non-relativistic (and non-Pauli) case from the points of view discussed here. 

We say that a system satisfies stability of the second kind or stability of matter 
if N^^En{Q) is bounded below independently of for all (open or in some cases 
sufficiently regular) i2 C M^. This is the version of stability mainly studied in |21 1. 

We will here use a slightly stronger notion which we refer to as grand canonical 
stability. We define the grand canonical ground state energy as 

E{|x,Q)=mf{S'{|x,\|/,A,E^_) y/e^, || v/|| = 1, A,Ex e Cq (R^M^^)}, 

where S'{jj,,y/,A,E±) was defined in (fTSI l. It of course depends on the species of 
particles. 

We say that a system is grand canonicaUy stable (with chemical potential jj.) if 

inf \n\-^E{ij.,n) > -oo. 

The infimum here is over all open sets £2 with bounded volume \Q \ (or possibly 
sufficiently regular sets if necessary, but we will not consider such cases here). 

The original proof of stability of matter is due to Dyson and Lenard |7 8 1 and 
later by a simpler method by Lieb and Thirring |28 1 . We will present a proof of grand 
canonical stability in a simple case relying on a combination of the two approaches. 

For grand canonically stable systems it is of interest to consider whether the 
thermodynamic limit 

lim \Q.\-^E{pL,Q.) (16) 

exists. The limit i2 — > K.^ can be given a precise meaning in different ways. Here we 
shall simply take the simple situation of the family of scaled copies LQ. of a fixed 
set Q. and let the real parameter L tend to infinity. 
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4.1 Stability of the first kind for non-relativistic particles 



We shall here prove the stability of the first kind for non-relativistic particles, i.e., 
particles with the kinetic energy (|4]i. 

Theorem 1 (Non-relativistic stability of the first kind). 

For all Iff e Cq (i2^)*'i-^" and all vector fields A, E G Cq (M^;R^) we have 

> -c\\v\\\ 

where the constant C > depends only on the number of particles N and their 
properties, i.e., on Vi , . . . , Vat S N, mi , . . . , > and Qi,..., Qn e M. 

This theorem follows easily from the diamagnetic inequahty and the Sobolev 
inequahty (see ||20| ). 

Theorem 2 (Diamagnetic Sobolev inequality). 

For all f e CJj[°(M^) and all A G C(~(R^;M^) there is a constant C > such that 

lj(-iv+m->livifif>c{lm'-)''\ 

An immediate corollary of this result (using simply Holder's inequality) is the fol- 
lowing bound on one-body Schrodinger operators. 

Corollary 1 (Lower bound on Schrodinger operator). 

For all f e Cq(E?), A e Cg (M^;M^), < Vi G L^/^(E?), and < y2 G L°°(R^) we 
have 

{f,{{-iV+Af-Vi-V2)f),2>-c[Jv^^' + \\V2\\S)\\f\\l2. 

We leave it to the reader to prove Theorem[T]from this corollary and the observation 
that |r|-l G L^/^(JS?) +L~(R3). 

The stability of the first kind holds even if the field energy 

is ignored. Moreover, it also holds if A is quantized, i.e., if we replace A{r) by the 
operator (fT3l i. This last statement follows since A (r) is a commuting family (indexed 
by r) and thus may be considered as a classical field. 
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4.2 Grand canonical stability 

We turn to the question of grand canonical stability. We will study this in the simple 
special case of two species of identical fermions with opposite charges. For grand 
canonical stability it is not necessary that all particles are fermions. It is, in fact, 
enough that all particles with one sign of the charge, i.e., say, all negatively charged 
particles form a collection of finitely many species of fermions. Stabihty of matter 
in this more general setting was proved in ||7][8]|28l (see also ITSl ) and the case of 
grand canonical stability was treated in 1 17|. 

One of the main ingredients in the proof of grand canonical stability is the use 
of the celebrated Lieb-Thirring inequality |28 | (see also 1211 ) which replaces the 
Sobolev inequaUty which we used in the proof of stability of the first kind. 

Theorem 3 (Lieb-Thirring inequality). 

Assume < V G L^I^{E?) and A e Co(K^;M^) then for all N we have on the anti- 
symmetric subspace — f\^[L?'{£l)Y the operator inequality 

t (^(-'V./+A(0))2-y(r,)) > -CmV2v J yV2, 

for a universal constant C > 0. In particular, it is independent of the number N of 
particles. 

Note the apparent similarity between the Lieb-Thirring inequality and the Corol- 
larylT] to the Sobolev inequality. The important difference is that Corollary[T]would 
only imply that 

t (i(-'V.+A(o))2-y(r,)) > -Cm^'N J V'l\ 

which, in fact, holds on all of (left as an exercise for the reader). The lower 
bound with a constant independent of holds only on the fermionic subspace. 

The Lieb-Thirring inequality relates the energy of a gas of independent parti- 
cles to the corresponding classical energy. The classical energy (ignoring internal 
degrees of freedom) would indeed be 

/ 1 ^{p^+A{r)f~V{r)drdp^-^-^my'-jv'l\ 

2L(p2+A(r))2-V(r)<0 

As a consequence of the Lieb-Thirring inequality we have the following lower 
bound on the kinetic energy of fermions confined to move in a bounded volume. 

Corollary 2. If the open set Q has finite volume \Q. \ then in /\^[L^(i2)]*' we have a 
universal constant C > such that 
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N 1 

L 5-(-«V^+A(r^))2 >Cm-^V-2/3AfV3|x2|-2/3. (17) 

Proof. If we use the Lieb-Thirring inequality with a constant potential V we obtain 

I ^ (-^'V,- +Airj)f >NV- Cn?l\Y'l^\Q \ 
which gives the estimate above after optimization in V. 

We now consider the situation with two species of identical non-relativistic fermions 
with masses m± > and charges ±2± where 2± > 0. For simplicity we assume that 
there are no internal degrees of freedom, i.e., v± = 1. In this case the Hamiltonian 
with particle numbers A'± for the two species is 

z.m+ j=N++l 
where the Coulomb energy is 

Q+Q- , ^ Ql . ^ Q- 



,7=li=A'+ + l I*"' 01 l<i<j<N+\''' 01 A'+<i'</<iV_+iV_ r' 01 

Note that we have numbered the positively charged particles l,...,N+ and the nega- 
tively charged particles A'^- + 1 , . . . , + . The Hamiltonian acts on the subspace 

/JV+ A'+ \ 

^= /\L2(i2)®/\L2(x2) nq°(i2^++^-). 



Theorem 4 (Simple case of grand canonical stability). 

The grand canonical energy in the finite volume set i2 C 

E{^,Q)=mf[{\ir,HN^^_\^)+H{N+N-) \ ^€9, \\w\\ = h 

£i,AeC2:(K^;K^)} 

satisfies the stability bound 

E{n,Q)>-C{n,m±,Q±)\Q\, 

for a constant C{n,m±,Q±) > depending only on H,m±,Q±. 

Proof. We define the distance from particle j to the nearest particle of the opposite 
charge, i.e.. 
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= Sj{ri,...,rN^+N-) 

_ J mini=N++i,...,N++N- \n -rj\, if j=l,...,N+ 

r^m=L....,N+ \n-rj\, if j =N+ + 1,,... ,N+ +N- ' 



Let Xj = ■^lB(ry,5^/2)' where B{rj,5j/2) denotes the ball centered at rj with radius 

8j/2 and lB{rj,Sj/2) is its characteristic function. Note that J Xj — ^■ 

We will use the following two observations: 
Observation 1: 

I I KT—n^L L Q+Q- ^;—^drdr 

i=\i=N+ + l \'' 'j\ i=\i=N+ + l •' •' y ' \ 



L T—Tl^ L 6+// ^^^^ 



Observation 2: 



1<1<./<A'+ I'"' Ol \<i<i<N^ 

and likewise for the -terms . 

The observations follow from Newton's Theorem: 



-1 



n?>^}v\<&n. "1 5-1(3 -4|r|25-2), if |r|< 5/2 



From the two observations above we arrive at the following lower bound on the 
Coulomb energy 

■p«p('-') ...... i2£e2^5-i_ll'^^^''- 

.7=1 ^ J=Af+ + l 

where we introduced the smeared charge density 



. W /■ PWp(0 , , , 12 ^2 .-1 12 ^2 .-1 



and used that 



;=i j=A'++i 



Using now the positive type (i.e., positivity of the Fourier transform) of the Coulomb 
kernel we find 

^ .7=1 )=A'+ + 1 

A similar application of the positive type of the Coulomb kernel goes back to an 
early paper of Onsager |[3Tll . who might have been the first to address the issue 
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of grand canonical stability. Better lower bounds on the Coulomb energy can be 
derived by more sophisticated use of the same ideas (see e.g. |in i29ll2Tl ). 
We are led to the following lower bound on the Hamiltonian 

Hn+ . n^ > Hn+ +Hn_+^i {N+ + N- ) 

where 

tiN. = E ^(-'^.- + Q+^^n)f - ^ I QW 

;=i + 7=1 

and likewise for . Observe now that for j — l,...,N+ the length 8j depends 
on the position rj and the positions r^^^i,. . . ,rN^+N_R of the negatively charged 
particles but not on the positions of the other positively charged particles. In other 
words we may write 

where 5{r) = min,=yv++i,...,A'++A'_ ki — ''I- We thus have a potential parameterized 
by the positions of the negatively charged particles. This observation allows us to 
use the Lieb-Thirring inequality Theorem [3] If we choose a parameter R (to be 
optimized over) and divide the space into the region where 5 (r) < R (a union of 
possibly intersecting balls of radius R) and 5{r) > R we obtain from the Lieb- 
Thirring inequality 



1 1 



2 



V J\r\<R 

> Cm;iA?+''^|X2 1-2/3 _ CQlf7iY^{N-R^/'- + \Q \R-^/^) 

where we saved half of the kinetic energy in the first inequality and estimated it by 
Corollary|2]in the second inequality. Finally, we optimized over the parameter R>0. 
Since the corresponding estimate holds for Hn_ we finally get the lower bound 

Hn^.N- > Cm;iA?+/^ 1X21-2/3+ Cm: 'A?!/^ 1X21-2/3 

-CQim'fN'J''\n\'/^ ^CQtm'J^N'flnl'/'' + ^{N++N-) 
> -C{l^,m±,Q±)\n\, 

where we have minimized in N±. We leave it to the reader to determine the exact 
form of the constant C{li,m±,Q±). 
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The same proof would work also if periodic external electric and magnetic fields 
were present, e.g., a situation describing a crystal structure. 
As should also be clear from the proof the field energy 




plays no role for stability in the present case. Moreover, as in the case discussed for 
stability of the first kind we could also have considered A quantized. 



4.3 Existence of the thermodynamic limit 

We will briefly discuss existence of the thermodynamic limit ( fT6] ). This was first 
proved by Lieb and Lebowitz 1221 for the case of several species of particles where 
all the species of, say, negatively charged particles are fermions. The method does 
not allow for an exterior periodic potential or magnetic field. In particular, the 
method does work in the case where the nuclei are confined to a periodic crystal 
arrangement. This case was later treated by Fefferman in |9|. In |,16,.17| an abstract 
method was developed to conclude existence of thermodynamic limits for Coulomb 
systems in great generality including periodic background potentials. 

Indeed, the method relies on establishing general abstract properties of the energy 
function that implies existence of the thermodynamic limit. 

We will just give a brief overview of the method. For the details and more precise 
definitions and assumptions we refer to |fT6llT7l . 

Let ^ ~ {£2 c M? open and bounded} and consider a map E : ^/M — > M with 
the following properties. Given a function a : [0,oo) with lim^^oo a(^) = 0, a subset 
M C ^ of sufficiently regular sets, constants fc, 5 > 0, and a reference set A e 
such that 

(Al) {^formalization). £(0) = 0. 

(A2) (Stability). \/Q £ E{Q) > ~K\n\. 

(A3) (Translation Invariance). \/Q ^ SP,,\/z ^I? , E{Q + z) = E{n). 

(A4) (Continuity). Vi2,i2' e ^, with Q.' <ZQ. and d{dQ.,dQ.') > 5, 
E{Q) <E{Q') + K\Q\Q'\ + \Q\a{\Q\). 

(A5) (Subaverage Property). For all Q. e we have 

E{n)>-^[ E{nng-{£A))dX{8)-\n\M() (i8) 

where dX is the Haar-measure of xi SO {3), (the group of isometrics of R^) and 
|i2|r := inf{|i2|, QCQ^ i2 e ^} is a regularized volume of i2 . 

If E : ^ — > M satisfies (A1-A5) then it is not very difficult to show that the 
thermodynamic limit 
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lim\£A\-^E{g£A) 

exists for all § e x: 5(9(3), i.e., it exists for all rotations or translations of the ref- 
erence set A. Under slightly more restrictive assumptions which we will not repeat 
here the limit holds for a very large class of regular sets. 

We see that (A2) is grand canonical stability. The difficult property to establish 
for Coulomb systems is (A5). For A being a simplex it is a consequence of the 
following result of Graf and Schenker lfT4l generalizing a somewhat similar estimate 
by Conlon, Lieb and Yau jS] : 

Theorem 5 (Graf-Schenker inequality). 

Let A be a simplex in M^. There exists a constant C such that for any N £ N, 
Qi,--,Qn e R ri,...,rAr e and any £ > 0, 

* 7=1 

This inequality follows by proving that the function 

Fir,r')^ [ l,(^ir,)l^eAirj)dX{g), 

Jr3>^SO(3) 

is of the form F(r, r') ^ g{\r — r'\) where g is such that kl^H^ ^^(kl) has positive 
Fourier transform. Recall that for a function / of positive type 

I G,e.//('-,-o)>-Ee//(o). 

l<i<j<N j=l 



5 Instability 

5.1 Examples of instability of the first kind 

As an example of a system that can show instability of the first kind we consider two 
relativistic particles with masses mi =m2 = l and charges Qi = — I and Q2 — Q>0. 
The kinetic energy is given by Q and we simply set A = 0. Thus the Hamiltonian 
is 

Ki -'"2I 

acting on the smooth compactly supported functions in L?(S?) L^{M?). Let 
ij/ e CJj[°(M'') be normalized, i.e., its square integral is one. and define ij/eirijfi) = 
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\j/{ri/£, nji) for ^ > 0. Note that i/^ is still normalized for all I. Then 



Thus if we let i tend to zero 



Q 



lim £{\if(,H\i/() = {wA 1 r W 

If Q is large enough we find that the right side is negative and hence for such a Q 

lini(v/£,//V^(') =-°o 

and the system is not stable of the first kind. 

On the other hand, if the negatively charged particles belong to a finite number 
of fermionic species and if the number of fermionic species, the maximal negative 
charge and the maximal positive charge satisfy appropriate bounds, then stability of 
matter holds l4l [121 12111^411291. 

A similar situation happens for non-relativistic particles interacting with mag- 
netic fields according to the Pauli operator Consider as an example a Hamilto- 
nian for two particles of mass m = 1 and charges Qi = Q> and Q2 = —Q <0- 

H{A) = Ua- (-/Vi - QA{n))f + Ua- (-/V2 + QA{r2))f 



2V V . - . 2' ' ^ ^ ' \n-r2\ 

1 



8;r 



|V®A| 



where we have chosen £ = (which is the energetically best choice). The instability 
in this case relies on the existence (see fTT, "301) of a non-zero \^ £ L^{M?) and a 
magnetic field A with / |V(g)Ap < 00 such that 

i(t7-HVi-A(n)))'r = o. 

We may assume that \j/ is normalized. If for £ > we set 

Vdri,r2)^i-^W{ri/e)w{r2/i) 

(which is also normalized) and Af{r) = {Q£)^^A{r/£) we obtain for the energy 
expectation 

\ \ri~r2\ I STtQ-^ J 

Again we see that if Q is large enough the right side is negative and hence for 
such a 2 we have as before hm^^„o {\l/e,H{Ai)y/() = —<^. As for the relativistic case 
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Stability of matter also holds in this case under appropriate conditions |fT0l[25 |. This 
problem with a quantized field has been treated in |l3]|TT|, the relativistic case with 
classical fields is considered in j26|, and the relativistic case with quantized field in 
i23J. 



5.2 Fermionic instability of the second kind 

As the final topic of these notes we will discuss instability of the second kind. 

We will first make a very simple general remark about instability of many-body 
systems with attractive interactions which has nothing to do with charged systems 
and holds even for fermions. 

Theorem 6 (Fermionic instability for attractive 2-body potentials). 

Assume that the potential W ; M" — > M satisfies W {r) < — c < for all r in a ball 
around the origin. Consider the N-body operator 

HN^t-Uj+ L W{r,-rj) 

j=l ^ l<i<j<N 

acting in the fermionic Hilbert space /\^L^(R"). If n > 3 then cannot be 
stable of the second kind, i.e., we can find a a sequence of normalized vectors 
y/N e L^{W) such that 

lim N^^ {xitnjHnWn) = 

Proof. Assume that < — c < on the ball of radius R centered at the origin. 
Define xj/^ as the (normalized) Slater determinant 

Wiru---,rN) = (W)-^/2det(«,(r,-))5=i 

where uj, j = 1,...,N are orthonormalized eigenfunctions corresponding to the 
lowest eigenvalues of the negative Laplacian with Dirichlet boundary conditions for 
the largest cube centered at the origin and contained in the ball of radius R. We 
extend the functions to be outside the cube. The functions uj are explicit and can 
be written in terms of sines and cosines. It is a simple straightforward calculation to 
show that in all dimensions n there is a constant C„ such that 

Comparing with Corollary |2] (written for the case n = 3) we see that there is always 
a similar lower bound. 
Thus 
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N-\xi/n,Hn\I/n) < CnN^I"R-^ - ^-{N- l)c. 
We see that instability occurs when n > 2. 
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5.3 Instability ofbosonic matter 



For matter consisting of charged particles we have discussed that the fermionic prop- 
erty ensures grand canonical stability. In this final section we will show that the 
fermionic property is indeed a necessity as stability fails for bosons. 

We consider two species of bosons with masses m± — 1, — —Q ~\, A = 
E± = 0. We describe them by the standard Schrodinger kinetic energy If we 
have positively charged particles and negatively charged particles we may 
write the Hamiltonian as 

N++N- 2 
j=l ^ l<i<j<N++N- I'' 'Jl 

where = 1 if ; = 1, . . . ,A?+ and ej = -1 if ; = + 1, . . . ,A^+ +A^_. The Hilbert 
space is 

M^.N^ =P+^(g)L2(R3)0P+_0L2(K3). 

This system is not stable of the second kind, in fact, the energy behaves like the 
number of particles to the 7/5-th power. The following precise asymptotics was 
conjectured by Dyson in [6 1. 

Theorem 7 (Dyson's formula). 

Let 

E{N) = inf mf{{w,HN^.N-W) I G n Co~(M^('^++'^-'), \M = U 

N++N-=N 

then as N ^ °° 



lim^^infji / |V0|2-/o/ 



0<<t>, I <t>^ = l}, (19) 



with Iq given by 

lo = (2/.)3/^/; 1 -.^ i.Uira. ^ ^^^^y (20) 

From the Sobolev inequality (Theorem|2]l it follows that the inf on the right of ( fT9] l 
is finite. In l&l Dyson proved an upper bound on E{N) of the form — cA^^/^ and thus 
indeed proved the instability of the second kind. In [5 1 a lower bound of the form 
—CN''/^ was established thus concluding that 7/5 is the correct power The theorem 
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was finally proved in ||27l [33l . In |fT9l Lieb proved that if the positively charged 
particles have infinite mass then the energy is much smaller, indeed, bounded above 
by —CN^/^ a corresponding lower bound had already been proved in ||7][8l- 

The proof of Theorem|7]relies on an application of Bogolubov's theory of super- 
fluidity llll. The charged system, in fact, forms a superfluid state. 

Dyson's formula ( fT9] l is proved by establishing the corresponding two inequal- 
ities. Establishing the lower bound is technically very involved and is beyond the 
scope of these notes. It is the content of the paper ll27l . We will here give a brief 
sketch of the proof of the upper bound from f33l. The upper bound is proved by 
finding an appropriate trial state. Here we are guided by Bogolubov's theory. 

It turns out that it is significantly easier to write down a grand canonical trial 
state than a canonical state. We are, however, interested in a canonical state. This 
will not be a serious problem as we will eventually be able to show that the state we 
construct is sharply peaked around the average particle number. We will ignore this 
point here and simply work with the grand canonical state. We refer the reader to 
Il33l for details. 

Another simplification is to consider the two species of bosons as one species 
with two internal degrees of freedom corresponding to the two signs of the charge. 
Constructing a trial state in this space will correspond to averaging over states with 
different numbers of positively and negatively charged particles. 

We are thus considering the Fock space ~ J?+(L^(R^)^). We write a func- 
tion / G L^(K^)^, as / = f{r,e), where e = ±1 is the sign of the charge. Let |0) be 
the vacuum vector in 

In constructing a bosonic trial state the first guess is to put all particles in the 
same one-particle state, i.e., to have a condensate. Let this state be represented by 
the (normalized) vector ^ e (R^ ) ^ . Introduce first the normaUzed grand canonical 
vector 



The corresponding state is an average over states with varying occupation in the 
condensate ^ . The average particle number in ^ is {E\a\{^ )a+ ) | S) = and the 
variance is also 



Thus this state is peaked around particle number with a standard deviation ^/N. 
There is a unitary operator U on ^+ such that 



Using this unitary we may also write |S) = f/|0). 

A pure condensate like this will however not give the correct state. It is impor- 
tant to build pair excitations too. This is achieved as follows. Let {/aj^^g 
orthonormal family in L^(R^)^ (they will represent the pair states). The normalized 




{Z\{a\{^)a+{^)f\Z)-{E\a\{^)a+{^)\E)^=N. 



U*a\{f)U^al{f) + VN{^Ja). 



Quantum Coulomb gases 27 
vector f e representing our trial state may be abstractly written 

= n (1 - ^a)'/'exp f E (<(/«) - ^^(^-/a))'') IS) (21) 

a=0 \a=0 ^ ^ ^ / 

= f/n(i-A^)i/4,^p/'£ _^,^(^^)2\ 

a=0 \a=0 ^ / 

We have introduced parameters < < 1 with Y^=q < to control the occu- 
pations in the pair states. For simplicity we will assume that ^ and {/aj^^o ^'"^ '^^'^ 
functions. 

We encode the information about the pair states in the positive semi-definite trace 
class operator on L^(M^)^ 

r=LT^I/«)(/«l- (22) 

a=0 ^ '^a 

In terms of this operator a lengthy but straightforward calculation shows that 

( -f^, « (/) - VA^(<^ , /) ) («+ (g) - VA^(g, ^ ) )f ) = (g, 7/> , 

(the inner product on the left is in ^+ and the one on the right is in L^(M^)^) and 
(f ,(fl^(/) - x/A^(^,/))'f') = 0. In particular-, 

('f^,<(/)«+(g)f) = (g,(A^|<^)(<^| + 7)/). 

Or equivalently using the field operators from Section [34l 

('f^,(/)(r,e)Xr',e')f)=A^<^(r,e)^(r',e') + 7('-,e;r',e') (23) 

where 7(r,e;r',e') is the integral kernel of 7. Likewise, 

^^^,a\{f)a\{g)^^) = (A^|^)(^h v/r(rTT)) /) , 

or 



('f^,(?)(r,e)>(r^^?0*'f')=A^^('-,e)'^('-^^0-v/l<r+T)('-,e;'-',e')- (24) 

Moreover, the state represented by f satisfies Wick's formula, which for the 4-point 
function reads 
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Here # refers to either a * (interpreted as complex conjugation on scalars) or no *. 
In particular, since ^ is real this gives 

f , {^{r,er - VN^(r,e)) (0(r',e')* - VN^{r',e')) 



= \Vriy+^)ir,e-y,e')\' + \r{r,e-y,e')\^ 

+Yir,e;r,e)Y{r\e';r\e') (25) 

Armed with these identities we can calculate the expectation of the energy in the 
state represented by f . 

First we will explain, for the special case of the charged Bose system, how to 
choose the condensate function ^ and the trace class operator 7. More precisely, we 
will specify their charge dependence. We set 

Ur,e) = ^Ur), (26) 

where ^0 is a real normalized function in L^{M? ). Thus the condensate function does 
not depend on the charge. The operator 7 on L? (M?)^ — {M? ) (g) will be chosen 
to have the form 

1/1 -1' 



r=7o^2V-i 1 

where 70 is a positive trace-class operator on L^(M^). Put differently, the integral 
kernel of 7 is chosen to be 

y{r,e;r\e')^^ee'yQ{r,r'). (27) 
The charge part of this operator is a rank one operator and thus we also have 

V7(7+l) = \/7o(7o + l)®^ (^_\ 

It is now straightforward to calculate the expectation of the Coulomb potential in 
the state represented by f . From ( l23ti27] i we obtain 



\ M=Ol<'<i<M 01 / 
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I [ [ee'Hr,er(p{r',e'r\r^r'\-'Hr,e)^{r',e')^ 



ee'=±l 



) 



A'Tr^2(K3) (^(70- V7o(7o + l))) • 



Here is the operator with integral kernel 



r,r') = ^oW|r-rr>^o(r'). 



The total energy expectation is 



( 



N+,N-=0 



>4/|v^.P 



+ iTr(-47o) +A^Tr(^(7o- V7o(7o + l))) 



The final step in the argument is to minimize the above expression over 70. More 
precisely this is done in a semiclassical approximation. We will only sketch this 
argument. The rigorous argument can again be found in lf33l . We assume that 70 is 
the quantization of a classical symbol f{r,p) > 0. The semiclassical approximation 
to the energy is then 



Minimizing this expression over /(r, p) and performing the p integration gives 



where Iq is given in ( l20l l. If we introduce the rescaling ^o('") = N^^^^'t'iN^^^r), 
where is also normalized then the energy expression above becomes 



which is exactly the expression conjectured by Dyson for the energy. 

Note that the instability is also reflected in the shrinking of the linear dimen- 
sion of the state with increasing A^. According to the scaling of ^0 above, the linear 
dimension of the state behaves like ~ N^^^^. 




r, p) + AnN\p I -2^0 {rf {fir, p) - ^ f{r,p){f{r,p) + \)) drdp. 
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